Quantum error correction of systematic errors using a quantum search framework 
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Composite pulses are a quantum control technique for canceling out systematic control errors. 
We present a new composite pulse sequence inspired by quantum search. Our technique can correct 
a wider variety of systematic errors - including, for example, nonlinear over-rotational errors - 
than previous techniques. Concatenation of the pulse sequence can reduce a systematic error to an 
arbitrarily small level. 
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I. INTRODUCTION 

Quantum error correction is perhaps the biggest hur- 
dle in building a quantum computer. Imperfect control 
operations are one of several sources of error. While error- 
correction schemes designed to correct for general errors 
Q no doubt also correct control errors, error-correction 
or error-avoidance schemes tuned to the dominant phys- 
ical error model are more efficient and practical. Spe- 
cialized error-correction schemes can also tolerate higher 
noise rates. 

This paper specializes to systematic control errors of 
the following form. When we try to apply the singlc- 
qubit pulse U — cxp(if?n ■ a), a 20 rotation about axis 
h, we in fact apply U = UV. The error is systematic 
in the sense that it is invertible; attempting to apply 
W in fact applies . The form of the error V is of 
course restricted. Previous authors 0,0 have considered 
the case of linear over-rotational errors: V = exp(ie(9n • 
(j) where e is fixed and small, but unknown. Here, we 
consider the case of general over-rotational errors, V = 
cxp[ie(|f?|, h)n ■ (?]. The amount of over-rotation, e(|0|,n), 
can now depend arbitrarily on the rotation angle 20 and 
also the axis of rotation h. 

Our error-correction method is a new composite pulse 
sequence, inspired by the generalization of quantum 
search known as amplitude amplification |J] . In this algo- 
rithm, the amplitude produced in a particular target sub- 
space by applying some unitary U to a source state is am- 
plified by successively repeating Ro(n)W Rt(n)U . Here 
-Ro(tt) and Rt(n) are selective reflections about the source 
and target, respectively. In standard quantum search, 
the source is |0"), U = H® n transverse Hadamard, and 
the target is a bit string |ar). In the subspace spanned by 
C7|O n ) and \x), the state vector steadily rotates toward 
| a;). Eventually, it rotates past the target. 

What happens in quantum search if we don't merely re- 



flect about the source and target, but instead add a phase 
other than 7r? It is well known that any phase bounded 
away from works to give a square-root speedup (with 
different constants); for example, this fact is used in a 
stronger form in Ambainis's element distinctness algo- 
rithm ||. One of us (LG) noticed that concatenating the 
basic sequence 

URo(ir/3)WR t (ir/3)U 

results in the state converging to the target subspace and 
not overshooting it, when viewed at times 3 fe , k € N 0|. 
Figure Ogives some geometrical intuition. 

For the present problem of systematic control errors, 
there is no source or target - we desire a "fully com- 
pensating" pulse sequence accurate on an arbitrary in- 
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FIG. 1: Looking down along the vector \s), the top di- 
agram shows why when {t\U\s} = (s|V|s) is close to 
one, {t\UR (7v/3)W Rt(Tr/3)U\s) = (s\VR V^ R V\s) is even 
closer to one. The bottom diagram gives similar intuition for 
a longer pulse sequence. 
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put - but a similar calculation still applies. We need 
merely choose a source arbitrarily, say the X +1 eigen- 
state |+), and set the target accordingly (to U\+)). As- 
sume n = (0,0,1). Let i?o = exp^X) a 7r/3 rotation 
about the X axis, and Rt = URqW . When we apply the 
sequence of noisy pulses 

RtURoWRtURl, 

the different systematic errors in both U and the noisy 
7r/3 rotations largely cancel out, leaving behind a higher- 
order error. (The extra pulses at cither end adjust for the 
phase difference between |+) and its orthogonal comple- 
ment |— ) which would otherwise be introduced.) This 
correction sequence can be concatenated on itself in a cer- 
tain way to reduce errors arbitrarily. A directly related 
method also applies to over-rotational errors in two-qubit 
gates .Therefore this composite pulse sequence allows 
for an arbitrarily accurate set of universal gates, giving 
a threshold result for this error model. 

We also consider another error model of systematic er- 
rors in even the rotation axis n: V = exp(ie-<r). Here the 
error e may depend on the rotation angle 28 but, except 
for a specific coordinate change, not on the axis h. For 
example, Rq and R t are related by the coordinate change 
Rt = URqW. We require that the errors be related by 
the same coordinate change, or R t = URqW. 

Section [H] describes the basic idea behind the compos- 
ite pulse sequence, by explaining its behavior when the 
7r/3 rotations are perfect. In the two following sections, 
we extend the error model to the two cases described 
above. 

Composite pulse sequences are an important, practical 
quantum control tool for removing systematic errors in a 
variety of quantum information processing implementa- 
tions 0, H .We need to show that our correction sequence 
remains practical. While the error models we address are 
more general than the linear over-rotational errors which 
have previously been considered, the control requirement 
is also stricter. We typically require the ability to rotate 
about an arbitrary axis in the Bloch sphere, not just 
one in the xy plane. If rotations arc only allowed about 
axes in the xy plane as in most NMR-type models, then 
our method only applies to correct 7r pulses. In Sec. [V] 
we compare our method, with it pulses and linear over- 
rotational errors, to previous fully-compensating compos- 
ite pulse sequences, particularly those recently discovered 



II. PERFECT tt/3 PULSES 

It is instructive to start with just the core idea of our 
composite pulse sequence, and build up the analysis from 
there. Consider the sequence 

U [x) = RlURoU^RtURl = URIvRoV^RoVRI, (1) 

where again V = U^U. Here we maintain a distinction 
between Rq = W RtU and Rq because the errors in the 



two terms might be different. For the rest of this section, 
however, assume Rq and Rt are perfect, so Rq = Rq = 
R = exp(ifX). 

Write V = exp(«e* • <?), where a = (X, Y, Z). Generally, 
each term of e may be nonzero. To measure the closeness 
of U^ x > to U, we compute a power series expansion of 
(Tr(A • tfuW),Tr(Y ■ WU^),Tr(Z ■ UW™)). Wc 
obtain 

(2*6, - iV3(e 2 y + e 2 z ) + 0(|ef ), 

2ie 3 z + 2ie z e 2 y + 0(\e\ 5 )). 

The first term is first order in e x because our correction 
Rq is a rotation about the x axis, and commutes with er- 
rors in the X direction. Errors in the Y and Z directions 
are symmetrically cancelled out, leaving only third-order 
terms. 

We can express this result quite simply. Assume e x ,y,z 
is an a, 6, cth order term. Then the X direction error 
order after the X direction composite pulse correction is 
applied is min{a, 26, 2c}, the Y direction error order is 
min{36, 6 + 2c} and symmetrically for the Z error. In 
shorthand, we write 

(a; 6; c) -> (a, 26, 2c; 36, 6 + 2c; 3c, c + 26). (2) 

The underset X here refers to X correction, and it is un- 
derstood that we take a minimum on each of the three 
terms on the right. This notation lets us quickly un- 
derstand what happens when we concatenate correction 
sequences. To concatenate when 7r/3 pulses are per- 
fect, just substitute the previous level's composite pulse 
sequence for U. A level k concatenation will require 
rife — 3rifc_i + 4 pulses, so the sequence length grows like 
4 fc . For example, starting with only first-order Z error, 
and applying an X correction gives 

(oo;oo; 1) (2;oo;3). 

At this point, it is best to apply a Y correction, since that 
cancels out errors in both X and Z directions. (Y cor- 
rection is symmetrical to X correction, except with 7r/3 
rotations about the y axis and the same axis conjugated 
by U.) At the next level of concatenation, Z correction 
will be optimal, and so on: 

(oo;oo;l) (2;co;3)^ (6; 4; 7) 

^ (14; 12; 7) (14; 26; 21) 

-> (42; 26; 49) (94; 78; 49). (3) 

After 3 6 = 729 pulses of U or W, and 1456 perfect tt/3 
pulses (about six axes), the error is only 0(|e| 49 ). 

Remark 1 (Generalization). The question of whether 
this pulse sequence generalizes deserves further study. 
We have been able to find a pulse sequence with five 
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applications of U or W, and perfect rotations by tt/5 
or 3tt/5, which on input |±) = ^(|0) ± |1)) achieves a 

fidelity error of 0(e w ): 

\(±\uWR»(%)WR t (-f)UR (-f)WR t (%)U\±)\ 2 

= l-0(e w ). (4) 

(See Fig. ^f 01 ' geometrical intuition.) However, this se- 
quence gives no improvement with imperfect correction 
rotations. 

Remark 2 (Error measurement). For us it is key 

to measure the direction of the error, as well as its 
magnitude. How does our method of measuring error 
compare to other reasonable methods? On a particu- 
lar input state, the difference in the fidelity from one is 
quadratically smaller than our measure. The so-called 
infidelity between U and U*, or 1 — i|Tr WfJ*\ is used in 
0- @ ,and is also quadratically smaller than our mea- 
sure. Brown ct al Q use as their measure of distance the 
trace distance Tr\U — U*\, which depends on the global 
phase of the operators. Our correction sequence does 
not give higher order accuracy in the global phase, but a 
simple modification of the trace distance optimizes over 
global phases, and then this measure of error is of the 
same order as ours. 



III. IMPERFECT tt/3 PULSES, ERROR 
ANGLE-DEPENDENT & AXIS-INDEPENDENT 



Let us consider the more realistic case that the tt/3 
pulses arc themselves erroneous. Assume that the error 
in a rotation depends only on the rotation angle, not on 
the rotation axis. When we attempt to apply exp(i9n-a), 
we actually apply exp{i[9 + e(9)]h ■ a}. Here e(9) can be 
an arbitrary function of 9, which is however always small 
[order 0(|e|)] . The error amount does not depend on the 
rotation axis n. Previous work has only considered the 
less-general case of linear errors, e(9) = e9. 

In fact, let us generalize our calculations slightly fur- 
ther. We will allow errors in the tt/3 pulses besides 
just over-rotation errors, except these errors must be the 
same under change of coordinates by U. In particular, 
Rt = URqW . It isn't clear in what physical models this 
will be an appropriate base error assumption - perhaps 
one in which the entire apparatus for applying a tt/3 ro- 
tation is rotated about the qubit in three dimensions as 
U acts in R 3 , or equivalently, the qubit is physically ro- 
tated. However, the added generality will be necessary 
for considering concatenation of this correction sequence. 

Write Rq = Rq exp(id ■ a). We obtain 



Tr(X • rffjW) = 2ie x + 2i(V3S y + S z )e y - iV3(e 2 y + e 2 ) - 2i(5 y - V3S z )e z + 0(|ef + |5||ef + |,5| 2 |el) 



Tr(r W) 



2i(V3S y + 6 z )e x - 4i(V3S x - V36 y S z + S 2 z )e y - V3i(V36 y + S z )e 2 y + 2ie^ 
-2i{V35 2 y - 26 y S z - V3S 2 z )e z - V3i(V36 y - S z )e 2 z + 2ie y e 2 z + 0(\e] 3 + |5| 2 |^), 



Tr(Z ■ Wlj( x ') can be determined by symmetry. In our shorthand notation, with 5 x ,y,z being d, e, /th order, respec- 
tively, 



a, e + b 7 f + b, 26, e + c, / + c, 2c; 
(a; b; c) -> I e + a, / + a, d + b, e + f + b, 2/ + b, e + 26, / + 2b, 36, 2e + c, e + / + c, 2/ + c, e + 2c, / + 2c, 6 + 2c; 

1 e + a, / + a, 2e + 6, e + / + 6, 2/ + 6, e + 26, / + 26, d + c, 2e + c, e + / + c, 26 + c, e + 2c, / + 2c, 3c 



r 



For example, taking d = e = f = oo, we recover Eq. (|2"|) 
from the perfect tt/3 pulse case. In the case of first-order 
over-rotation, d = 1, e = / = oo, 



(oo;oo;l) x -^(4;4;4). 



(0) 



To obtain arbitrarily accurate rotations, it is most effec- 
tive to correct both the applications of U and the tt/3 
correction pulses. So at this point correct the tt/3 pulses 
until d = e = f = 4. Note that applying such a cor- 
rection maintains the invariant that the error in a pulse 
depend only on the angle and not the axis. Now 



(4; 4; 4) — > (12; 12; 12). 



(7) 



Every three levels of concatenation (both on U and the 
tt/3 pulses) increases the error order by a factor of three. 
Therefore obtaining error tolerance to a desired amount 
e* requires poly- logarithmically many pulses in 1/e*. 



IV. IMPERFECT tt/3 PULSES, ERROR BOTH 
ANGLE- AND AXIS-DEPENDENT 

What if the error in tt/3 pulses depends on which basis 
they are carried out in, i.e., R t ^ URqU^. Can wc still 
obtain arbitrarily accurate pulses? Perhaps surprisingly, 
the answer is yes, if the error is of a restricted form: 
only over-rotation errors. However, the orders will not 
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grow exponentially quickly in the number of concatena- 
tion levels, only linearly, implying that error tolerance to 
an amount will require polynomially many pulses in 
1/e* instead of only poly-logarithmically many. 

Write R = Roexp(i5X), R t = UR exp{iSX)W. Ex- 
panding (Tr(X-WU^),Tr(Y-WU^),Tr(Z-WU^)), 
we obtain 

/ 2iex-iV3(e 2 y +e 2 z )+0(\8\e 2 + \e] 3 ), \ 

V -2i^(8+S)£ z -2i(8-S)e v +2 t 4+2ie z t 2 y +0(S 2 \e\ + \i] 5 ) J 

assuming |<5| = 0(|<5|). In our shorthand notation, with 
5 first order, the rule is 



(a; 6; c) 



x 



(a, 26, 2c; 36, 6 + 2c, 1 + 6, 1 + c; 3c, c + 26, 1 + 6, 1 + c). 



For example, 

(oo; oo; 1) 



^ (2; 2; 2) ^(2; 3; 3) 
^ (3;3;3)^(3;4;4), 



(8) 



(9) 



and so on, with every two levels of concatenation increas- 
ing the error order by one. Note that we do not concate- 
nate corrections onto the 7r/3 pulses, because then the 
error would no longer be simply over-rotational. (Even 
with a more general expansion, it turns out that the con- 
vergence is still only be linear in the number of concate- 
nation levels.) 



V. 7T PULSES IN NMR 

While our method corrects against more general types 
of errors than previous composite pulse sequences, it also 
has a stronger requirement. Namely, we must be able to 
apply a 7r/3 rotation about the x, y and z axes, and also 
about those same axes in the [/-transformed basis. In 
most current proposed quantum information implemen- 
tations, primitive rotations are only allowed about axes 
in the xy plane. For U, Rt and Rq all to be rotations 
about axes in the xy plane, it must be that U is a rota- 
tion by an integer multiple of ir. This is a considerable 
restriction on the applicability of our method. Still, for 7r 
pulses our correction succeeds in a setting more general 
than that for which previous methods could correct; for 
example, we can correct for non- linear over-rotations. 

Assume now that the systematic error is in fact a linear 
over-rotation; when we try to apply cxjp(i8n - a), we actu- 
ally apply cxp[i#(l+e)n-er]. Assume U — exp(i : |A) = X . 
Our 7r/3 correction method, correcting in the Y direction, 
leaves behind second-order errors, with a sequence length 
of 37r + 47r/3. (Practical composite pulse sequences need 
to be as short as possible, in order to minimize any non- 
ideal effects not accounted for in our error model.) We 
cannot concatenate a Z correction, but can concatenate 
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FIG. 2: Infidelity plotted against the base linear over-rotation 
parameter e. Solid line: unprotected tt pulse. Dashed line: 
B2 and B4 pulse sequences. Dotted line: n/3 Y correction on 
the 7r pulse, and concatenated onto a symmetrized B2 (BB1) 
pulse sequence. 



alternately X and Y corrections to reduce the error arbi- 
trarily. Concatenating an X correction onto the Y cor- 
rection leaves behind a third-order error, with a sequence 
length of (14i)7r. 

How does the 7r/3 correction method compare with 
previous correction methods? The most practical correc- 
tion methods previously known were the B2 (also known 
as BB1) correction sequence, and the recently discovered 
B4 sequence Q- These sequences are implemented as 
follows: 



B2: 



B4: 



^ 1 (7r) J R 3 ^(27r)^ 1 ( 7 r))C/ (10) 

(R^(ir)R^ 2 (2Tr)R^(n)) 4 \ 
,(-27r)iL*,(-47r)i2fc(-27i-)) U, (11) 

x(^ 2 (7r)^30 2 (27r)^ 2 (7T))V 



where R^O) = exp(z e ^ ^ (cos (j)X + sin</>F)), cos^i = 

— j^, and cos ^2 = — 2~f^- They leave behind third- and 
fifth- order errors, respectively. The total sequence length 
for correcting a 7r rotation is 57r for B2 and 4l7r for 
B4. Hence our method, in this linear over-rotation error 
model, seems to offer little over the plain B2 sequence. 
Table|Ilgives values for the infidelities 1- ||Tr UW^ | for 
various correction sequences. (As previously remarked, 
the infidelity is quadratically smaller than the trace dis- 
tance, so the infidelity error orders for the three possibil- 
ities 7r/3 Y, B2 and B4 are 4, 6 and 10, respectively.) 

Docs our method complement previous correction 
methods? To answer this question, we must determine 
the direction of the error left behind after a correction se- 
quence. For example, B2 and B4 are each implemented as 
BiU, where Bi is some particular pulse sequence not in- 
volving U. A simple calculation shows that both B2 and 
B4 leave behind an error which is has relatively large X 
and Y components. Therefore, concatenating on X or Y 
correction will not increase the error order. We can how- 
ever find an axis in the xy plane which is approximately 
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orthogonal to the xy component of the error, and correct 
along this axis. Alternatively, we can symmetrize the B2 
and B4 sequences into exp[ij(l+e)X]Bi exp[ij(l+e)X]. 
In this more symmetrical form, the error magnitude is un- 
changed, but the direction is entirely into the xz plane. 
Therefore, we can simply apply a Y correction to the 
symmetrized sequences. Table [i] compares the infidelities 
of Y correction concatenated onto the symmetrized B2 
and B4 correction sequences. Note that the former case 
gives fourth-order protection with a sequence length of 
only (16|)7r; this gives a new, perhaps practical, com- 
promise between the B2 and B4 correction sequences. 
Figure |2 plots the fidelities for e > 0. 

VI. CONCLUSION 

We have presented two main results. The 7r/3 cor- 
rection sequence protects against general errors which 
depend arbitrarily on the rotation angle but not the 



rotation axis. The same sequence protects against over- 
rotational error which depends arbitrarily on both the 
rotation axis and the angle of rotation. Previously, com- 
posite pulse correction sequences were only known for 
the cases when the error was independent of the rotation 
axis, and depended linearly on the rotation angle. 

Moreover, our composite pulse correction sequence 
concatenates nicely to reduce errors arbitrarily. In 
the first case, the overhead number of pulses is poly- 
logarithmic in the desired accuracy, and in the second 
case the overhead is polynomial. 

However, our correction sequence in general requires 
primitive rotations about arbitrary axes in the Bloch 
sphere, and only applies to correct it pulses in the typ- 
ical situation in which rotations are only allowed about 
axes in the xy plane. For correcting it pulses, our method 
concatenated on top of a B2 pulse correction provides a 
new compromise between B2 and B4. 

B. R. acknowledges support from NSF ITR Grant 
CCR-0121555, and ARO Grant DAAD 19-03-1-0082. 
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TABLE I: Infidelities 1 - ±|Tr U^U W \ of naive and variously compensated exp(i^X) pulses. 
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